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Abstract 

In a certain class of differential-difference equations for dissipative systems, we 
show that hyperbolic tangent model is the only the nonlinear system of equations 
which can admit some particular solutions of the Toda lattice. We give one pa- 
rameter family of exact solutions, which include as special cases the Toda lattice 
solutions as well as the Whitham's solutions in the NewelPs model. Our solutions 
can be used to describe temporal-spatial density patterns observed in the optimal 
velocity model for traffic flow. 
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1 Introduction 



This paper concerns analytic description of temporal-spatial density patterns generated 
in a class of discrete nonlinear systems of dynamical equations: 

G[x n (t), x n (t), x n (t + r), x n (t + r)] = V[Ax n (t)}, (n = 1, 2, N), (1.1) 

where the r.h.s. depends only on the forward difference between the discrete indices Ax n (t) = 
x n -i(t) —x n (t). It would be interesting if we could find exact solutions for such dissipative 
systems. 

It is the purpose of this paper to show that there is actually a nonlinear system of 
differential-difference equations for which we may find exact solutions. To show that, 
we shall construct nonlinear equations from given candidates for exact solutions instead 
of trying to solve a given system of G[- ■ •] and V(Ax) exactly. This is analogous with 
Toda's way of finding his lattice model[|l] 0). A clue which led him to the discovery^ 
is the observation that a specific combination of logarithm of theta functions obeys a 
simple system of equations of motion for an exponentially interacting lattice of particles. 
The particular solutions, referred to as Toda solutions hereafter, correspond to a periodic 
wave |l] as well as a solitary wave[@]. These are used here to describe temporal-spatial 
density patterns, and to construct nonlinear systems of equations. 

Requiring V(Ax) to be single- valued, we show that only the first-order differential- 
difference equations |§||. 

x n {t + r) = V[Ax n (t)}, (1.2) 
with a hyperbolic tangent (tanh) OV function V(Ax) can admit the Toda solutions^ (See 

(ill))- 

The first-order equations ( |1.2j ) and its truncated form of the second-order differential 
equations with "friction terms" |J 1| 

rx n (t) + x n (t) = V[Az n (t)]. (1.3) 

have been extensively discussed in the optimal velocity (OV) model||][0|||7]] to describe 
the formation of density pattern in a congested flow of traffic; they may be also relevant 
to the density patterns in granular materials^. To be concrete, we employ the OV model 
terminology: the basic coordinate x n and the forward difference Ax n correspond to the 
position of nth car and its headway, respectively. The latter is the distance between the 
car and the preceding (n — l)th car. The idea postulating eqs. (|1.2| ) is that a driver adjusts 
the car's velocity x n (t) according to the observed headway Ax n . The delay time r is the 
time lag for the driver and car to reach the optimal velocity V(Ax) when the traffic flow 
is changed. In the OV model, spontaneously appears a density pattern with regions of 

^ee ref.g for review articles. 

2 Although the functional form of the Toda solutions here is the same as the original one in the Toda 
lattice, the dispersion relations of the two systems are obviously different. We use the term of "Toda 
solutions" up to the dispersion relations. 

3 See, for example, ref.B- 
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high density where cars move slowly, and low-density regions where the velocity of cars 
is high. When these two regions occur alternatively in a lane, a spatial-temporal pattern 
is generated, which may be viewed as a density wave. What we show here is that the 
pattern can be described by the Toda solutions of the nonlinear equations ([2.151) . 

We also show that the our tanh model can admit one parameter family of exact 
solutions which contains a width parameter S. The Toda solution for a "cnoidal wave" 
consists of sharp pulses with a fixed width <5 T . In contrast, our solutions describe pulses 
with arbitrary width. Unless 5 = 5t, the solutions do not obey the Toda equation The 
solution for S = 15t with an integer I may be interpreted to describe a "bound state" of 
solitons. It is certainly intriguing that such solutions exist for dissipative systems. We 
have confirmed stability of the solutions via a computer simulation. 

The nonlinear systemF] of eqs.( |1.2j ) was first discussed by Newell [f|]. He gave exact 
solutions for r = in the model with Vn(Ax) — &i[l — exp{— 6 2 (Ax — 63)}]. Some years 
ago, Whitham]^] obtained Toda-like exact solutions of the Newell's model for r^O. He 
found that a relation between parameters which specify a wave propagation is crucial 
for the existence of the exact solutions. In our construction, this relation, which we call 
the Whitham relation, is necessary for the OV function to be single-valued. We will 
see that the Whitham's solutions belong to our family of solutions: they have a width 
5 = 5w = 1/2 — &r/2, and our tanh OV function reduces to the Newell's function in this 
case. 

This paper is organized as follows. The next section describes construction of the 
model equations by using the Toda solution for a "cnoidal wave" . In section 3, we present 
an ansatz for one parameter family of exact solutions as an extension of the Toda solutions, 
and discuss the stability of the exact solutions. We obtain a solitary wave solution in the 
large modulus limit in section 4. Summary and a few remarks on the exact solutions are 
given in the final section. Appendix lists some formulae used in this paper. 



2 Construction of the tanh model 
2.1 Toda solution for periodic wave 

In a simulation of the OV model, we may generate a stable pattern of a congested flow. 
There, the pulses with almost the same maximal values are observed in the velocity as 
well as the headway of a car. Here we consider the density pattern analytically by using 
the Toda solution corresponding to a "cnoidal wave"]]]]. For the position of the nth car, 
it is given by 



^0 


(yt- 


-s) 


00 ( 


vt- 


n+l\ 
A ) 



x n (t) = A In — ^ + Ct -nh, (2.1) 



4 Some analytic solutions were found for the system fll.3| ) with piecewise linear OV functions |9| |n| , or 
in an asymptotic method to investigate the long time behavior in the vicinity of a critical point. 
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where A, X, v, C, h are constants. #o(^) is the theta function (See ( |A.4|) in Appendix). 
The headway is then given by 



Ax n (t) = A 



In 



tfo (id - f ) (id - a±l 



+ /i. 



(2.2) 



To write equations given below in compact forms, we introduce the notations: 



v 

So 



n 



Vo = x 



sn(2Kvo), cq = cn(2Kvo), do = dn(2Kvo), 



(2.3) 



where sn, cn and dn are the Jacobian elliptic functions, K = K(k) is the complete 
elliptic integrals of the first kind. All the Jacobian elliptic functions used in this paper 
are assumed to have a modulus k. 



Using formula (|A.1| ) in Appendix which relates the theta functions with the elliptic 
functions, one obtains 



where 



X 



k 2 s 2 sn 2 {2Kv) 
Ax n -p 



1 - e 



x 



A 



(3 = h + 2A In 



^o(^o) 

MO) ' 



(2.4) 



(2.5) 



Therefore, the elliptic functions, sn(2Kv), cn(2Kv) and dn(2Kv), are expressed in terms 
of e x . The velocity and acceleration of the nth car are given by 



x n {t) = 2AKv[Z{2Kv) - Z {2K(v - v ))\ + C 
x n {t) = AA(Ku) 2 \dn 2 (2Kv) - dn 2 (2K(v - v )) 



[2.6) 



where Z(u) is the Jacobian zeta function. The formulae ( |A.2| ), ( |A.6| ) and ( |A.7| ) in Ap- 
pendix lead to 



x n (t) = 7 - (2AKuk 2 so)sn(2Kv) 
7 = 2AKv Z(2Kv ) + C, 



c d ■ sn(2Kv) - s ■ cn(2Kv)dn(2Kv) 
l-k 2 s%-sn 2 (2Kv) 



and 

x n (t) = AA(Ku) 2 



dn 2 (2Kv) - 



do ■ dn(2Kv) + k 2 s c ■ sn(2Kv)cn(2Kv)V 



l-k 2 s 2 -sn 2 (2Kv) 
The velocity and acceleration at a time t + r are calculated to be 



(2.7) 



(2.8) 



x n (t + r) = 2AKu[Z(2K(v + v 1 )) - Z(2K(v + v 2 ))} + C 
x n (t + r) = 4A(Kv) 2 \da 2 (2K(v + v 1 )) - dn 2 (2K(v + v 2 )) 



(2.9) 
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where V\ = ut, v 2 = v\ — Vq. These are again written by the elliptic functions, 

"cidi • sn(2Kv) + si ■ cn(2Kv)dn{2Kv) 



± n (t + r) = 7 - (2AKuk 2 s ) 



and 



x n (t + r) = AA(Kv 



1 - k 2 s{ ■ sn 2 (2Kv) 
c 2 d 2 ■ sn(2Kv) + s 2 ■ cn(2Kv)dn(2Kv)~ 
1 - k 2 s\ ■ sn 2 {2Kv) J ' 

dt ■ dn(2Kv) - A; 2 Sl Ci • sn(2Kv)cn(2Kv^ ^ 2 



1 - k 2 s\ ■ sn 2 (2Kv) 
d 2 ■ dn(2Kv) - k 2 s 2 c 2 ■ sn(2Kv)cn(2Kv^ ^ - 
1 - k 2 s\ ■ sn 2 (2Kv) 



where 



(2-10) 



Si = sn(2Kvi), Cj = cn(2 Kvi), di = dn(2Kvi) (i = 1,2). 
Eqs.(|2.7|), ( |2.8| ), ( |2.10| ), ( |2.11| ) are used to construct the dynamical equations. 



(2.11) 



(2.12) 



2.2 Uniqueness of the model 

Here under appropriate assumptions we show that only the x n (t + r) is the fundamental 
element to be incorporated into the dynamical equations. 

Let us rename the velocity and acceleration variables as, 
(2/1)2/2)1/3)2/4) = (£ n (t), x n (t), x n (t + r), x n (t + r)). In general, y's are all double-valued 
functions of Ax n because of the presence of terms proportional to sn(2Kv) -cn(2Kv). One 
may impose some cancellation conditions which make their coefficients to vanish. We first 
examine if it is possible to construct single- valued dynamical equations without imposing 
such conditions. One might find a polynomial F = YUfi {Ui) ni in f° ur variables, where /'s 
are chosen in such a way that F becomes a single- valued function of Ax n , F = V(Ax n ). 
It gives a system of single- valued equations, but takes of the form: F(yi, y 2 , y 3 , y±; Ax n ) = 
V(Ax n ), where the l.h.s. necessarily depends on Ax n via the non- vanishing coefficients, 
/'s. The class of these models which have the Toda solution may not be out of our 
interest, but does not belong to the system ( |l.l| )p|. 

From the above results we are tempted to turn to the possibility of having single- valued 
y l by imposing cancellation conditions on the parameters, (s m , c m , d m ) (m = 0, 1, 2). They 
are functions of (i/t, A) and k. As discussed in the next section, the periodic boundary 
condition xn+i = X\ requires that A is order of the total number of the discrete elements 
(cars in the OV model), N. Therefore, it is not unreasonable to assume, say A > 2. This 
excludes Cq = 0. It is also obvious that Sq ^ 0, otherwise Ax n (t) = constant. In the same 
manner it can be shown that the cancellation conditions necessary for x n (t + r) to be 



5 This observation may apply to non-polynomial combinations of ?/' which are single-valued functions 
of Ax n . 
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single- valued are not acceptable. Thus, x n (t), x n (t), x n (t + r) cannot be single- valued, 
and are excluded from entering into the desired equations. 

One is then only left with the cancellation condition to make x n (t + r) single- valued^: 

s 1 c 2 d 2 + s 2 cidi = 0. (2.13) 

Solution for the condition, t> 2 = —vi, is essentially unique and acceptable. It leads to the 
Whitham relation |J 

2u\t = I. (2.14) 

One finds that x n (t + r) subject to the Whitham relation becomes only the variable which 
may be used to construct dynamical equations. For the first-order differential-difference 
equations, one obtains 



x n (t + r) = V(Ax n ) 

■Ax n - p 



£ + r^tanh 



2A 



(2.15) 



where £, rj and p are some constants determined by the parameters in the Toda solution 
fl2.ip . The tanh OV function has been extensively used in computer simulations [0] for 
the second-order differential equations x n = a[V(Ax n )—x n ]. In the traffic flow application, 
we have shown here the system of first-order differential-difference equations with the tanh 
OV function is only the model which admits the Toda solution. 

The following remarks are in order. 
(l)The nonlinear systems of equations as 

x n (t) = a(V[Ax n (t)]-x n (t)), (2.16) 
x n (t + r) = a(V[Ax n (t)]-x n (t)), (2.17) 

cannot admit the Toda solutions for single-valued OV functions. In other words, the OV 
functions become necessarily double-valued in order for the systems ( |2.16|) to have the 



Toda solutions. Each OV function may take two different values for a fixed headway, 
depending whether the car is accelerating or decelerating. Although it would not be 
unrealistic to have the double-valued optimal velocities, the simplicity of the OV model is 
lost. As pointed out in the above discussion, a way out of the appearance of the double- 
valued functions is to consider the equations like F{%)i, Ax n ) = V[Ax n }. For example, if the 
parameter a in ( J2.16p is allowed to have an appropriate Ax n dependence, the OV function 



V in the modified system of equations may become single-valued. Another possibility is to 
introduce the backward difference Ax n+ \ as well as the forward difference Ax n . The OV 
function becomes then a single-valued function of two variables. We shall not consider 



here such a class of the models|12|, though it will be intriguing. 



(2) The Whitham relation ( 2.14]) is found to be crucial for the existence of the exact 



6 For single- valued x n (t + r), x n (t) may be a double- valued function of Ax n upon the formation of a 
congested flow. 
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solutions. The relation is interpreted as follows. In a density wave propagation, the 
headway may be written in the form 

Ax n (t) = f(vt-j)=f(v) 

Arc n _ 1 (t) = f(v(t + T)-j^=f(v + vT), T = (2.18) 

which implies that the nth headway has the same time dependence as that of the {n — l)th 
one apart from a time lag r ||10|| . The Whitham relation ( [2.1 4|) then takes the form 

T = 2 r. (2.19) 

Imposing this relation, one obtains a differential-difference equation for a single universal 
function f{v) 

y±f(v) = V[f(v + u(T-r))} - V [f(v - vr)} 

= V[f{v + vr)] - V[f{v-ur)}. (2.20) 

Note that the argument of the left in fl2.20| ) is midway between those on the right. It 
suggests that a lattice model interpretation may be possible for the present system once 
a temporal-spatial pattern is generated. 



3 Ansatz for one parameter family of solutions 



The "cnoidal wave" described by the Toda solution (|2.1| ) consists of periodic pulses. There 
is a parameter, A, which adjusts the period of the pulses. In the Toda solution, however, 
the width of each pulse is fixed. In this section, we present one parameter family of 
solutions which consist of pulses with arbitrary width. 

In the context of the OV model, the motivation for constructing such solutions is 
explained as follows: In a computer simulation, the density pattern is generated under 
suitable conditions. It is described by the alternative appearance of the low-density (high- 
velocity) regions and the high-density (low- velocity) regions in a traffic flow. A congested 
region contains a bunch of cars. For the Toda solution, take a circuit of length L, and 
impose the periodic boundary condition, xn+i = x\- Then, the parameter A is fixed to 
be A = N/m, where m is the number of the congested regions in the circuit. The solution 
( |2.1|) with A > (A < 0) describes a car which spends most of its time in congested (free) 
regions. Thus, the distance between a pair of a kink and an anti-kink in x n in the Toda 
solution is much shorter than that observed in the simulation; the pair looks like a sharp 
pulse. It is desirable therefore to explore analytic solutions which consist of "trapezoidal 
pulses" with larger width. 

Our ansatz for one parameter family of solutions is given by 



x n {t)=A In — ) f + Ct -nh, (3.1) 
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where 5 is related to the size of a congested or free region. Without loss of generality, 
the range of 5 is restricted as < 5 < 1/20. One may observe that 5 = 5t = 1/(2A) 
corresponds to the Toda solution. Let us emphasize that the ansatz otherwise does not 
solve the Toda equation 

Let us see that this ansatz does solve the tanh model ( |2.15| ). The ansatz ( |3.1| ) slightly 
modifies the expression for the headway, 



e x -l 



k 2 sn 2 (2Kv) 



e x ■ sn 2 2K(5 - i) - sn 2 2K(5 + i 



(3.2) 



where 



X 



AXr 



A 







Let us introduce the following notations: 



h + 2A In 



1 



Ms - ±) 



(3.3) 



P = sn\2K5) - sn 2 (2K(5- — )), 

2A 

Q = -sn 2 (2K6) + sn 2 (2K(5 + —)), 

2X 

s s = sn(2K5), c s = cn{2K5), d s = dn{2KS) 
For r which satisfies the relation ( p.!4| ), 

k 2 sa 2 {2Kv) 



x n (t + r) = AAKv 



Z{2K5) — sscsds 



+ C. 



(3.4) 



(3.5) 



1 - k 2 s 2 s ■ sn 2 (2Kv) 

Taking the same step as before, we obtain the following OV function for the present case, 



V(Ax) = AAKv Z(2K5) + C + (AAKv ■ s 5 c s d 5 ) 



P-e x + Q 



(3.6) 



It may be rewritten in a more familiar form: 



V(Ax) = C + 2AKv 



2Z{2K5) + s s c s d. 



5 • 



P Q 



+ 2AKv s s c s d s + 1 



/ Ax-J3 1 Q" 

tanh in — 

\ 2A 2 P 



(3.7) 



This completes the proof that the ansatz ( |3.1| ) solves ( |2.15| ). 

Note that the OV function is described by a tanh function only when P/Q > or 
in the range, 1/(4A) < 5 < 1/2 — 1/(4A), outside of which it becomes a coth function; 

7 One may take 5 to be — 1/2 < S < 1/2 because of the periodicity of the ^o-function. Furthermore, 
the difference between x n for S = —Sq and that for 5 = 1/2 — 6q can be absorbed into a constant shift of 
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on the boundaries it is exponential functions. In particular, it should be noted that if 
one takes 5 = 5w = 1/2 — 1/(4A) so that Q = 0, the solution ( |3.1|) reproduces the 
Whitham's solution ||, and the OV function in ( |3.7| ) reduces to the one in the Newell 
model, V(Ax) = 6^1 - exp{-6 2 (Ax - 63)}]. 

Comparing the OV function in fl2.15| ) with that of ( |3.7| ), we have 



C + 2AKv 



2Z(2KS) + s s c s d s 



P 



Q. 



2AKv s s c s d 5 ■ f i + , 

$o(5-l/(2\)) + Aln p- 



(3.1 



It is crucially important to see the stability of our "trapezoidal pulse" solutions. We 
have confirmed it via a numerical study. 

Suppose iV cars run on the lane with length L according to the equations of motion 
given with the time delay r and the OV function. The linear stability analysis tells us 
that above some critical value for r, we find a range of density, or the average headway 
h = L/N, for which a homogeneous flow is linearly unstable. What we have learned 
from numerical studies are: (1) once a congested flow is fully developed, a car repeats 
the behaviour of its preceding one with the time delay T and the entire pattern moves 
backward with some fixed velocity fBfllQfl; (2) the trajectory of a car on the headway- 
velocity plane crosses the OV function at two points well outside the unstable region and 
the line connecting the two points has the slope ~ 1/T and passes through the inflection 
point. Obviously, the T and vb are the outputs of the simulation. 

In our analytical study reported in this paper, the set of parameters (N, L,t,£, p, rf) as 
well as (T,vb) for the global pattern are written with the parameters in our "trapezoidal 
pulse" solutions. This implies that there are nontrivial relations among the quantities 
stated above. Among others, the relation, T = 2r, cannot be foreseen from numerical 
study and should be understood as a prediction of our analytical study. 

We have performed a simulation to compare with our analytical results. To start a 
simulation we need to choose appropriate values for r, N and L. Since they are written 
with parameters in our "trapezoidal pulse" solutions, one needs to search a set of pa- 
rameters in the solutions which give natural values for the three quantities as those for 
the OV model. The following conditions are important: (1) the value for r allows the 
linearly unstable region; (2) h is in the unstable regionP]. The simulation with the chosen 
values for r, N and L gave us a hysteresis loop on the headway-velocity plane. This is 
to be compared with the loop obtained from our analytical solutions.0 As clearly shown 

8 Since we started from the solutions, logically we do not have any strong reason for the presence of 
appropriate set of parameters which allows an interpretation in the context of the OV model. 

9 Here are parameters for the simulation: r = 0.582, N = A = 10 , h = 1.89; the OV function is 
V(Ax) = tanh 2 + tanh(Ax - 2). The result is to be compared with eq.(U) for 2K 5 = 3, k 2 = 0.99999 
and C = 0.866. The DODAM package is used for solving differential equations. 
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in Fig.l, they are actually found to be the same! In particular, the prediction from the 
solution, T = 2r, is confirmed. 



velocity 



l 




5 10 15 20 



Fig.l: 

Velocities from a simulation and "trapezoidal pulse" solution; 
dots are from the simulation and the line denotes the solution. 

Therefore we conclude that our "trapezoidal pulse" solutions are the exact stable 
solutions for the OV model, which describe the fully developed congested flows. 



4 A solitary wave solution 

We have discussed so far the periodic solution, and constructed the model with the cor- 
responding exact solution. In the Toda lattice, it has been known that a solitary wave 
solution [g] is obtained from the periodic wave solution Q by taking the large modulus 
limit k — > 1. We use the same procedure here. In the solution (|3.1| ), the exponential of 
Ax n is related to dn 2 (2Kv ). It consists of an infinite number of pulses, each of which 
is expressed by square of the hyperbolic secant (See the decomposition formula (|A.8|) in 
Appendix). A solitary wave limit is defined by k — > 1 (K — > oo), v —>■ , A — > oo and 

5 —>■ oo with 2Ku = ur, 2K/ A = 1/Xr and 2K5 = 5r held fixed. In this limit, only one of 
the pulses survives to describe a solitary wave, while all the others disappear. Actually, 
since 

hm ^° {{VR + 6r ~ Vm /^I 2K ^ - C ° Sh ( VR + 6r ~ V ~^) (41) 
fc-i $ ((vr :- 5 R - v 0R /2)/2K) cos h L R _ § R ' ' 

where v R = u R t - n/X R and v 0R = 1/X R , x n is given by 

cosh (v R + 5 R - V -^) 

L — - 

cosh (v R - 5r - 



2 



cosh (v R + 6 R - V -^) 
x n (t) = A In ) U- + Ct ~nh. (4.2) 
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Therefore, x n (t) consists of a pair of a kink and an anti-kink; x n (t) has two pulses with 
different sign: 



x n (t) 
x n (t) 



tanh [v R - 5 R 



Av 



R 



_ vqr \ 
2 J 

sech 2 (vr + 5 r - ^ ) - sech 2 (v R -5 R - ] 



ta.nl) j Vr + Sr- ^Y^j 



+ a 



2 J 



2 J 



(4.3) 



This implies that the kink and the anti-kink is separated by the distance 25 r. One can 
directly confirm that ( [4.2| ) and ( (4.3| ) provide another particular solution of ( |2.15| ). 



5 Summary and outlook 

We have given a method to construct the dynamical equations for many-body dissipative 
systems, which admit the Toda lattice solutions describing periodic as well as solitary 
density waves. The uniqueness of the tanh model of the first-order differential- difference 
equations is shown. It should be emphasized again that the Whitham relation plays an 
important role. It is remarkable that the model admits one parameter family of exact 
solutions which are not present in the exact solutions of the Toda lattice. 

For the width parameter 5 = 5\y — 1/2 — 1/(4 A), our solutions become the Whitham's 
solutions, and the tanh OV function reduces to the Newell's function. As described above, 
we do know there are ranges for parameters to make our solutions stable. However it is not 
obvious to us for the moment the parameters for the Whitham's solutions fall into these 
ranges. So further study via computer simulation as analytical method will be needed to 
clarify this point. 

The model considered here seems to be not exactly solvable, yet carries some nature 
to form patterns stably. This might suggest that there are important classes of "partly 
solvable models" , sharing stability of solutions with solitons in the exactly solvable models. 
They may serve to reveal the underlying universal nature of pattern formations. In this 
connection, it is a highly non-trivial task to investigate soliton picture in the present 
model. In the Toda lattice, solitons are stable and independent except the time interval 
during collisions, and multi-soliton states are essentially described by a succession of two- 
soliton collisions. Soliton picture associated with our model seems to be different from 
this. What we have observed in simulations is as follows: when two congested regions 
collide, they make a single large congested region; after a certain time, some number 
of such large congested regions move with the same velocity and the entire pattern on 
the lane does not change. So we expect that solitons in this model are stable against 
perturbations, but may not necessarily "solitary": suppose that two solitons run in the 
same direction and one of them catches up the other, they would form one-soliton state, 
their "bound state", whose size is given by a sum of those for two solitons. Our new 
solutions discussed in the previous section could be understood as "bound states" . Given 
such a picture, it would be very important and interesting to know how two solitons 
would behave in this model. In order to make the soliton picture clear, analysis of the 
two-soliton sector would therefore be crucial. 
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From our results, it is tempting to speculate that there are dissipative systems which 
acquire (at least, partial) integrability when temporal-spatial patterns are formed: there 
could be some crucial relations similar to the Whitham's for other system to recover such 
properties. 
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A Appendix 

Here we list useful formulae^ used in the text . 

iM^L - In «^-21n^W=ln[l-^sn=(2^)sn^(2^)], (A.l) 



Z(-u) — 7i(w) = Z(u — w) — k 2 snu snw sn(u — w), (A. 2) 

ni \ m \ crv, \ 2 fc 2 snw cnw dnw sn 2 M 

Z(u + w) - Z(u-w)-2Z(w) = — , (A.3) 

1 — k z sn z w sn z u 

where K is the complete elliptic integral of the first kind. The theta function do(v) and 
the Jacobian zeta function Z(u) are defined by 

oo oo 

Mv) = Qo 11(1- 2g 2n - 1 cos2vrt; + 2 An ~ 2 ), q = - q 2n ) (A.4) 

n=l n=l 

Z(u) = ^o(^), Z'(ti) = dn\ - |, (A.5) 

with q = e~ 7rK '/ K . K' is the complete elliptic integral of the first kind for the compli- 
mentary modulus k' = yl — k 2 , and E the complete elliptic integral of the second kind. 
Other formulae are: 

snu cnw dnw; — snw cnu dnu 

m[u-w) = j^- , (A.6) 

1 — k z sn z u sn z w 

dnu dnw + k 2 snu cnu snw emu 
dn(u-w) = — — , A.7 

1 — fc^sru sn z w 



and 

2 oo 



dnH2Kv) - | = (JL) ± -h 2 \^(v-l) 



l=—oo 



IT 



2KK' 
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See ref. || for formulae of the elliptic functions. 



11 



References 

[1] M. Toda, J. Phys. Soc. Japan 22 431 (1967). 
[2] M. Toda, J. Phys. Soc. Japan 23 501 (1967). 

[3] M. Toda, Prog. Theor. Phys. Suppl. 45 174 (1970); Phys. Rep. 18 1 (1975). 
[4] G.F. Newell, Oper. Res. 9 209 (1961). 

[5] G.B. Whitham, Proc. Royal Soc. London, A 428, 49 (1990) 

[6] M. Bando, K. Hasebe, A. Nakayama, A. Shibata and Y. Sugiyama, Phys. Rev E51, 
1035 (1995). 

[7] M. Bando, K. Hasebe, A. Nakayama, A. Shibata and Y. Sugiyama, Jpn. J. Ind. Appl. 
Math. 11, 203 (1994). 

[8] B. S. Kerner and P. Konhauser, Phys. Rev. E48 (1993) 2335; Phys. Rev. E50 (1994) 
54; 

S. Sasa and H. Hayakawa, Euor. Phys. Lett. 17, 685 (1992) 
[9] Y. Sugiyama and H. Yamada, Phys. Rev. E 55 (1997) 7749. 
[10] K. Nakanishi, K. Itoh, Y. Igarashi and M. Bando, Phys. Rev E55, 6519 (1997). 
[11] K. Komatsu and S. Sasa, Phys. Rev E52, 5574 (1995). 
[12] H. Hayakawa and K. Nakanishi, Phys. Rev E57, 3839 (1998). 



12 



